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We re-analyze the observed magnitude-redshift relation of type la supernovae (SNe la) 
and examine the possibility that the apparent acceleration of the cosmic expansion is not 
caused by dark energy but is instead a consequence of the large-scale inhomogeneities in the 
universe. We propose a method to phenomenologically describe the effects of the large-scale 
inhomogeneities without relying on the specific toy models of the inhomogeneous universe. 
This method clearly illustrates how the post-Friedmannian effects of inhomogeneities, i.e. 
the effects due to the deviation from a perfectly homogeneous and isotropic model, act as 
an effective cosmological constant in the magnitude-redshift relation of SNe la. 
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The Cosmological Principle, which states that our universe is described by the 
homogeneous and isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) model, 
in some averaged sense, is a working hypothesis which has been widely accepted 
among cosmologists. The present, past, and future evolution of the FLRW model is 
determined by a few constant parameters, such as the Hubble parameter, Ho, the 
matter density parameter, f2 m , and the cosmological constant, A (or the normalized 
parameter Q/y = A/3Hq). The determination of the cosmological parameters is one 
of the main purposes of observational cosmology. 

The recent observations of type la supernovae (SNe la) now strongly suggest 
the acceleration of the cosmic expansion^®'®'® As long as we employ perfectly 
$_i ■ homogeneous and isotropic FLRW models, this requires dark energy, an exotic energy 



component which accelerates the cosmic expansion with its negative pressure. 

Instead of introducing such a mysterious energy component, there have been at- 
tempts to explain the apparent accelerated expansion of the universe resulting from 
the large-scale inhomogeneities in the universe. For example, TomitaP®'^®® us- 
ing his local void model, and Iguchi et al.j^ using the Lemaitre-Tolman-BondPS'^'^ 
(LTB) model, studied the possibility of explaining the observed magnitude-redshift 
(m-z) relation of SNe la. Moreover, recently Alnes et alP^ have concluded that not 
only the m-z relation of SNe la but also the position of the first peak in the cosmic 
microwave background (CMB) anisotropy can be explained by the inhomogeneity in 
the LTB model. However, these works depend specifically on simplified toy models. 
Therefore, due to the lack of strong support for such toy models as providing realistic 
descriptions of our universe, the study of inhomogeneous effects in the universe is 
not the mainstream of research in cosmology. 
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In this article, we re-analyze the observed m-z relation of SNe la and point 
out some theoretical possibilities of the inhomogeneity interpretation to explain the 
apparent acceleration. Then, we propose a method to phenomenologically describe 
the effects of the large-scale inhomogeneities in the universe, without relying on 
specific toy models. This method clearly illustrates how the post-Friedmannian 
effects of inhomogeneities, i.e. the effects due to the deviation from a perfectly 
homogeneous and isotropic FLRW model, act effectively as a cosmological constant 
in the magnitude-redshift relation of SNe la. 

§2. The magnitude-redshift relation of SNe la 

The apparent magnitude m of a SN la of absolute magnitude M, at redshift z, 

is 

where Dl(z) is the luminosity distance in units of parsecs. The luminosity distance 
is obtained by solving the propagation of light ray bundles through space-time. In 
the FLRW universe, it is written in the forirP^'^ 

n < \ c(l + z) 

d l{z) 



Hq^/I — f2 m — fi A 



dz' 



x sinh Jl-Q m - Q A / = . (2-2) 

V -A) \/{l + V m z')(l + z' f -z'(2 + z>) Q A J 

The luminosity distance D^{z) is a slightly complicated function of z with three 
constant parameters, Hq, f2 m , and fi A . 

As an illustration, we use the observed SNe la data presented in the paper of 
Perlmutter et alP In total, 60 SNe la with redshift in the range 0.014 < z < 0.830 
are listed in Tables 1 and 2 of Ref. [2]). Because all of the SNe la have redshifts 
satisfying z < 1, the luminosity distance Dl(z) may be most usefully expressed as a 
power series, 

D L (z) = ^(z + d 2 z 2 + d 3 z 3 + ■■■), (2-3) 
where the expansion coefficients d 2 and d 3 are given by^ 1 

d 2 = j(2-n m + 2n A ), (2-4) 

d 3 = I {Mm 2 + 4 n A 2 - 4 Q m Q A - 2 Q m - 4 (2 A ) . (2-5) 



Substituting Eq. (|2 3j) into the m-z relation Eq. (|2-ip gives 

m = M-5 + 51og 10 D L (z) 

= M + 51og 10 (z + d 2 z 2 + d 3 z 3 ) , (2-6) 

where M. = M — 5 + 51og 10 c/ 'Hq. This quantity is often called as the "Hubble- 
constant-free absolute magnitude'® or the "magnitude zero-point" Note that 
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Taylor expansion up to at least 0(z 3 ) is necessary to determine the three parameters 
Hq, f2 m , and Qa by fitting to the data. Once the best fit coefficients d 2 and 0I3 are 
obtained from such a fitting, we can calculate the cosmological parameters fl m and 
SI A as follows: 

rt m = 2(1 -da) (1-2 (fe) -2 da, (2-7) 
f2 A = d 2 (2d 2 -l)-d 3 . (2-8) 




et al.P with the parameter values fi m = 0.28 and Qa = 0.72. 

In Fig. 1, we plot the Hubble diagram for all 60 SNe la. Also plotted there is 
the best fit m-z curve for a flat cosmology with the parameter values Q m = 0.28 and 
SI A = 0.72, which was determined by Perlmutter et al.^ 

§3. Inhomogeneous interpretation? 

In order to formulate an alternative interpretation of the m-z relation without 
dark energy, it is instructive to divide the whole SNe la data set into two parts. For 
the sake of convenience, we define SNe la with redshifts satisfying z < 0.2 as low- 
redshift (low-z) SNe la, and those with redshifts satisfying z > 0.3 as high-redshift 
(high-z) SNe la. Then, the low-z data set consists of 20 SNe la with redshifts in the 
range 0.014 < z < 0.18, and the high-z data set consists of 40 SNe la with redshifts 
in the range 0.320 < z < 0.830. No SNe la are in the range 0.2 < z < 0.3. 
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redshift z 

Fig. 2. Hubble diagram for the 20 low-redshift (z < 0.2) SNe la data set and the best fit ra-z curve 
of a zero-yl cosmology with the parameter values M = 24.01 and J? m = 0.26. 
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redshift z 

Fig. 3. Hubble diagram for the 40 high-redshift (z > 0.3) SNe la data set and the best fit m-z 
curve for a zero-yl cosmology with the parameter values M — 24.27 and f2 m = 0.35. 
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In this paper, we use an implementation of the nonlinear least-square methods 
of the Levenberg-Marquardt algorithm^ :Q2} to fit the m-z relation Eq. (|2-6p to each 
data set. 

Figure 2 displays only the low- z SNe la data set, along with the best fit m-z 
curve of a zero-yl cosmology with the parameter values Ai = 24.01 and fl m = 0.26. 
Since the cosmological constant A does not play an important role for z *C 1, it is 
evident that the zero-yl cosmology can fit the low-z data set fairly well. 

Figure 3 displays only the high-z SNe la data set, along with the best fit m-z 
curve of a zero-/l cosmology with the parameter values Ai = 24.27 and fl m = 0.35. 
It is interesting that a zero-yl cosmology can also accurately fit the high-z SNe la 
data set, although the values of the cosmological parameters Ai and fi m in that case 
are different from the best fit values for the \ow-z SNe la. It should also be noted 
that the high-z best fit m-z curve (Ai = 24.27, Q m = 0.35, Qa = 0) behaves very 
similarly to that of the flat cosmology with Q m = 0.28 and Qa = 0.72 obtained by 
Perlmutter et alPJ 1 in the range 0.3 < z < 1. 




redshift z 

Fig. 4. Hubble diagram for the 20 low-redshift (z < 0.2) SNe la and the 40 high-redshift (z > 0.3) 
SNe la. Also plotted are the best fit m-z curves of a zero-/l cosmology for each data set and the 
best fit flat cosmology with Q m — 0.28 and Qa = 0.72 for the entire SNe la data set obtained 
by Perlmutter et alP as a reference. 

Figure 4 summarizes all of the above results. There it is seen that an open 
zero-yl FLRW model with the parameter values Ai = 24.01 and f2 m = 0.26 can fit 
the low-z SNe la data, whereas another open zero-yl FLRW model with Ai = 24.27 
and Q m = 0.35 can fit the high-z data. No cosmological constant nor dark energy is 
necessary to fit either part. However, a positive cosmological constant {Qa = 0.72) 



6 



M. Kasai 



is necessary to fit the entire SNe la data set with a single flat FLRW cosmology. It is 
evident from Fig. 4 that this conclusion simply comes from the asymptotic behavior 
of the theoretical curve with fl m = 0.28 and Qa = 0.72 and does not depend on the 
particular choices of the SNe la data sets. 

Because M. = M — 5 + 51og 10 c/Hq, where the Hubble distance c/Hq is in units 
of parsecs, the result that different values of M. fit different redshift data sets implies 
the following possibilities: 

1. The absolute magnitude M of the high-z SNe la is systematically different from 
that of the low-z SNe la. 

2. The speed of light c is different in different redshift regions. 

3. Hq in the high-z region is slightly different from that in the low-z region. 

Let us examine the third possibility, namely, the inhomogeneity interpretation. 
If the difference between M{iow-z) = 24.01 and .M(high-z) = 24.27 is due to the 
inhomogeneity of the Hubble parameter Hq, we can estimate the following ratio 
representing the inhomogeneity in Hq between the two different redshift regions as 

24.01 - 24.27 = 51og 10 °\ 8 . (3-1) 

This yields Ho(high-z) = 0.89 Hq(\ow-z). An inhomogeneity in which the value of Hq 
in the high-z region is 11% smaller than that in the low- 2; region may be sufficient 
to explain the observed m-z relation for SNe la, without the need to introduce a 
cosmological constant or dark energy. 

§4. Effects of large-scale inhomogeneities on the luminosity distance 

In the previous section, we pointed out the interesting possibility that it may be 
possible to account for the observed m-z relation for SNe la by the large-scale in- 
homogeneities in the universe, without introducing a cosmological constant or dark 
energy. In order to examine this possibility further, we need to study the prop- 
agation of light ray bundles through the inhomogeneous universe and obtain the 
luminosity distance Dl{z) as a function of the redshift z. Since the actual space- 
time inhomogeneities of the present universe are not known in detail, the usual 
approach employs simplified toy models of the inhomogeneous universe. For exam- 
ple, TomitsP'^-EJ®'® used a local void model, and Iguchi et alP* used the LTB 
model. 

We take another approach. In this section, we propose a method to phenomeno- 
logically describe the effects of the large-scale inhomogeneities on the luminosity dis- 
tance, without relying on specific toy models of the inhomogeneous universe. This 
method clearly illustrates how the "post-Friedmannian" effects of inhomogeneities, 
i.e. the effects due to the deviation from a perfectly homogeneous and isotropic 
FLRW model, act effectively as a cosmological constant in the magnitude-redshift 
relation of SNe la. 

In the perfectly homogeneous FLRW models, Hq denotes the expansion rate at 
the present time, to, and it is constant over the entire t = to hypersurface, due to the 
perfect spatial homogeneity. In inhomogeneous universes, however, the expansion 
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rate is naturally dependent on the spatial positions. Therefore, Hq is not constant 
and may depend on z: 

H H (z). (4-1) 
Analogously, the density parameter Q m may also depend on z: 

n m => n m {z). (4-2) 

In general inhomogeneous universes, these cosmological parameters may also 
depend on the angular direction due to spatial anisotropics. Already in 1966, Kris- 
tian and Sachs^J emphasized the importance of observing angular variations in the 
various cosmological effects. Kasai and SasakP^ and KasaP^ derived formulae for 
cosmological observations in a linearly perturbed FLRW model in gauge-invariant 
manner and found the existence of a quadrupole anisotropy of the Hubble parame- 
ter Hq, which is directly proportional to the gauge-invariant scalar potential. They 
found that "the perturbed space-time behaves as a Friedmann-like universe with 
the direction-dependent Hq and qo"-^ I n this paper, however, we concentrate on 
the z dependence of the cosmological observables. Future investigations will include 
consideration of the angular dependences of Hq and Q m . 

In the region z < 1, the cosmological parameters can be expressed in the power 
series forms 

Hq{z) = Ho(l + h 1 z + h 2 z 2 + ■■■) , (4-3) 
f2 m (z) = Q m (1 + wiz + uj 2 z 2 H ), (4-4) 

where Hq = Hq{z = 0) and Q m = f2 m (z = 0), and the expansion coefficients 
hi,h,2i ... , uji,ll>2, ■ ■ ■ represent the "post-Friedmannian" corrections due to spatial 
inhomogeneities. The models reduce to the FLRW if and only if all of these co- 
efficients vanish. Under the assumption that Hq(z) and f} m {z) are slowly varying 
functions of z, we substitute them for Hq and Q m in Eq. (|2 2[) and obtain the fol- 
lowing power series formula: 

D L (z) = 4~( z + ^2 z 2 + d 3 z 3 + --- ), (4-5) 
-no 

where the expansion coefficients are now 

d 2 = - (2 - ti m + 2 n A ) - hi, (4-6) 



da = i (pi, + 4 Qa 2 - 4 fl m Q A - 2 ti m - 4 Q A } 

+1 U{h x ) 2 - ±hi ~ ±h 2 - hi {2Q A -Q m -2)- Q m \ . (4-7) 

Now we can illustrate how the "post-Friedmannian" corrections of the spatial 
inhomogeneities act effectively as a cosmological constant. Suppose that astronomers 
obtain the best fit parameters d 2 and d% by fitting the luminosity distance formula 
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Eq. (|4-5p to the observed m-z data of SNe la. If they assume that the universe is 
homogeneous and isotropic, they will simply use Eqs. (|2-7p and (|2-8|) to calculate the 
cosmological parameters. Then, even if the true value of the cosmological constant 
is zero, they will obtain the following effective value for the cosmological constant 
from Eq. fitty : 

Qf = d 2 (2d 2 -l) -d 3 

= X - U>h x n m + |wiAn -2/ii+4 (h x f + 4fr 2 J (4-8) 

This clearly shows that the "post-Friedmannian" correction terms hi,h 2 , and 0J\ to- 
gether act effectively as a cosmological constant. The effective value fi c ^ is unrelated 
to the true value of the cosmological constant. It simply results from the erroneous 
assumption that the universe is perfectly homogeneous and that Hq and fl m are 
constant on the t = to hypersurface. 

In the same way, we also obtain the following effective value for the density 
parameter from Eq. ()2-7|) : 

= 2(l- d 2 ) (l - 2 d 2 ) - 2 d 3 

= (l + h x + iwij An + 3/ii + 2(/n) 2 + 2/i 2 . (4-9) 

Following the procedure described in §2., we obtain the best fit values for the 
cosmological parameters as (2^ = 0.28 and = 0.72. Unfortunately, however, 
these do not completely determine the inhomogeneities of the actual universe. This 
is simply because the two cosmological parameters and J2^ ff are functions of four 
parameters D m ,hi,h 2 , and u)\ (or possibly five, including QjCj. The data fitting of 
the m-z relation itself only yields constraints on some sets of the post-Friedmannian 
parameters, but it does not determine completely the values of each parameter in- 
dependently. For this reason, it is highly desirable to incorporate other independent 
observations, such as CMB data, gravitational lensing data, and so on, in order to 
determine the extent to which our universe is homogeneous or inhomogeneous. 

§5. Summary 

We have re-analyzed the observed m-z relation of SNe la proposed by Perlmutter 
et alP and have examined the possibility that the apparent acceleration of the 
cosmic expansion is a consequence of large-scale inhomogeneities in the universe. As 
previously found by Perlmutter et al.P a positive cosmological constant is necessary 
to fit the whole data set, consisting of 60 SNe la in the redshift range 0.014 < z < 
0.830, with a single FLRW model. They obtained the best fit values Q m = 0.28 and 
17/1 = 0.72 for a flat cosmology. 

In order to examine the feasibility of the inhomogeneity interpretation, we di- 
vided the SNe la data into two parts, low-z and high- z data sets. The low- z (z < 0.2) 
data set consists of 20 SNe la in the redshift range 0.014 < z < 0.18, and the high- 
z (z > 0.3) data set consists of 40 SNe la in the range 0.320 < z < 0.830. We 
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were able to fit the low-z and high-z data sets, respectively, with two zero- A FLRW 
cosmologies. The best fit parameters are M = 24.01 and Q m = 0.26 for the low-z 
data set and A4 = 24.27 fi m = 0.35 for the high-z data set. The difference between 
the values of M implies that the Hubble parameter Hq in the high-z: region is 11% 
smaller than that in the low-z region. This indicates the possibility of the cosmolog- 
ical parameters Ho and Q m being dependent on the redshift z. It also suggests that 
the nearby low-z region is a local void, i.e. a less dense and more rapidly expanding 
than the outer high-z region. This local void interpretation is consistent with the 
results of previous works employing inhomogeneous toy models EKE^E}© MM ,03} 

Inspired by the above results from the data fittings, we proposed a method to 
phenomenologically describe the effects of the large-scale inhomogeneities, without 
relying on specific toy models of the inhomogeneous universe. In general inhomoge- 
neous universes, the Hubble expansion rate Hq and the density parameter f} m are 
naturally dependent on the spatial coordinates, and therefore the redshift z. We 
expanded H (z) and f2 m (z) into power series with respect to z and obtained a lumi- 
nosity distance formula with corrections due to the large-scale inhomogeneities. This 
distance formula clearly illustrates how the corrections resulting from the large-scale 
inhomogeneities act effectively as a cosmological constant. 

In this paper, we re-analyzed only the SNe la data presented by Perlmutter et 
alP Other data sets, including those with z > 1 SNe la, should also be examined. 
It is also noted that the m-z relation of SNe la itself does not completely determine 
the inhomogeneities of the actual universe. The data fitting of the m-z relation itself 
only gives constraints on some sets of the post-Friedmannian parameters. It does not 
completely determine the values of each parameter independently. For this reason, 
it is highly desirable to incorporate other independent observations, such as CMB 
data, gravitational lensing data, and so on, in order to determine the extent to which 
our universe is homogeneous or inhomogeneous. 
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